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An  Efficient  Technique  for  the  Calculation 
of  Velocity  —  Acceleration  Periodograms 


The  velocity  —  acceleration  periodograra  associated  with  the  (complex) 

data  samples  r  ,  ...  r„  ,  is  defined  by 
o  M-l 


H*,a) 


N-l 

-x 

k=0 


J2*fkA  j2na(kA)£ 


(1) 


where  A  denotes  the  (uniform)  time  separation  between  successive  data  points. 

-1  -2 

P  is  periodic  in  f  with  period  A  and  periodic  in  a  with  period  A  so  that 

P  need  only  be  evaluated  over  the  (f,a)  region  defined  by  0  <  f  <  A-^, 

_2 

0  <  a  <  A  .  Furthermore,  since  the  velocity  and  acceleration  resolutions  of 

-1  -2 

the  periodogram  are  givm (approximately)  by  (NA)  and  (NA)  respectively, 
it  is  usually  sufficient  to  evaluate  P  at  the  discrete  points  given  by 
f  =  n(NA)~^,  a  =  m(NA)~2  where  n  =  0,  1,  ...  N-l  and  m  =  0,  ...  N^-l.  These 
considerations  transform  the  original  periodogram  problem  to  the  evaluation 
of  the  expression: 


N-l  2 

F(n,m)  =  £  V1”* 

k=0 


(2) 


where  n  =  0,  1,  ...  N-l;  m  =  0,  1, 


. ..  N2-!;  W  =  exp( j  2*/N),  V  =  exp(j  2«/N2). 
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Following  Cooley  and  Tukey*,  we  assume  that  N  =  2?  and  proceed  to  express 


the  integers  k,  n,  m  in  binary  form  as  follows: 


k  =  k  12p"1  +  ...  +  le.2  +  k 
p-1  1  c 


n  =  n  -2P^+...+n,  2  +  n 
p-1  1  o 


m  =  nig  ^22p  ^  +  ...  +  m^2  +  in. 


where  k^,  n^  and  take  on  the  values  0  and  1.  In  addition,  it  will  be 

2 

convenient  to  express  k  in  the  form 


k2  =  (k2)  +  ...  +  (k2) 

'p-1  '  'o 


2  2 

where  (k  ^  =  those  terms  in  k  that  depend  on  k^  ^  but  not  on  k 
k  , .  Thus, 

p-1 


(k  )  .  =  k 

p-iJ  p 


.  sf**1  f  k  a1”8  +  k2  ,  22^ 

-&  p-g  P-* 


g=^+l 


The  derivation  of  this  last  formula  is  straightforward  exercise.  Note  that 


*  Cooley  and  Tukey,  An  Algorithm  for  the  Machine  Calculation  of  Complex 

Fourier  Series,  Math,  of  Comp.  !£;  April,  1965. 
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except  when  l  =  p. 


(10^  ^  contains  a  factor 


Next  we  note  that 


W1* 


w 


(n  +  ...  n  ,  2p-1)(k  + 

o  p-1  y  x  o 


,.  +  k  ,  a5"1) 

p-i 


k  .  2p  ^n  .  0p-^/  _■> 

=  W  P-1  °  W%-2  ^  (no  +  nl2) 


-2/ 


W 


k  (n  +  . 
o  o 


'•  +  Vi 2P_1) 


and 


v“(kC 


2)  _  v^ra0  +  ••• +  ^p-i 2  p  1)t(k  )Q  +  ••• +  (k  )p_i^ 


(k^)  ,  (m  +  ...  +  m  ,  2P_1)  ... 

x  /p-lx  o  p— 1 


=  v 


v(k2)1(mo  +  ...  +  nt,^  22p'3)  v(k2)o(mo  +  ...  +  m^  2^“x) 


,2p-l> 


because  the  exponent  of  W  need  only  be  computed  modulo  N  =  2P  and  the  exponent 

2  2d 

of  V  need  only  be  computer  modulo  N  =  2  . 

With  some  obvious  changes  of  notation,  equation  (2)  now  can  be  written  in 
the  form 


1  w 
P-1 


m  , 
o* 
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•  •  • 


“I 


w 


k  (n  + 
o  o 


n  i  a5-1)  (k^)  ( 

P-1  „  O' 


/m.  +  ...  +  nu  2^ 
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,  2p"1n  (k^)  (m  +  ...  +m  .  2P“1) 

-1  o  p-lN  O  P“1 


-1, 


r<v  • 


Vi>  w 


P-1 


(3) 


For  computational  purposes,  it  is  convenient  to  think  of  equation  (3)  as  a 
sequence  of  p  calculations  as  follows:  First  compute 


w  * 


1c  « «  n  •  m  •  • 
p-2*  cr  cr 


•  •  m  i  ) 

P-1 


$  l**  - 


k  -  ^  (k^)  -  (m  +  ., 

k  ,)  wP-1  0  V  P'1  ° 

P-1 


.  m 


(k> 


p-i 


then  successively  compute  from  P^  Z  =  Z,  ...  p-1,  according  to  the 
formula 


P/V  **•  kp-i-l,  V  •••  nz- 1>  V  **•  mp+i-2^  ~ 


* 
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^  ^-l^V  '*•  V 

k  . 

P-i 


...  n&_2,  raQ,  ...  rap+x_3) 


W 


P-* 


-,p -!>t  ^^-l^  t -  2> 


(nQ  +  ...  ni_1  2  )  (k  )  ^  (m0  +  •” 
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Finally,  P  is  computed  from  the  formula, 
P 


+  V-f-2  ^ 


+£-2 


(5) 


P  (n  ,  ., 
p  o 


n  ,  m  ,  . . .  nu  n  ) 
p-1  o'  ^p-1' 


-i 


P  ..  (k  ,  n 
p-1  o  c 


n  o>  m  ) 
p-2'  o' 


%-3> 


o(no  +  •"  Vl  2P  l)  „(k‘Vmo  +  •••  +  “20-1 
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2p-l> 
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(6) 


The  last  computed  function  P  is  the  desired  function  P  given  by  equation  (3)* 

1? 

A  straightforward  computation  of  the  periodogram  from  equation  (2)  would 

O 

require  (N-l)lr  computations.  (A  computation  is  defined  as  being  the  perform¬ 
ance  of  two  complex  multiplications  followed  by  an  addition.  Thus,  each 
evaluation  of  the  sum  in  equation  (2)  requires  N-l  computations  and,  since 
there  are  N • values  of  n  and  m  for  which  the  sum  must  be  evaluated,  the 

O 

resulting  number  of  computations  is  (H-l)N  • )  The  computation  method  just 
proposed  requires  many  fewer  computations  as  will  now  be  demonstrated. 

The  calculation  of  P^  requires  2*2^  =  2^  computations  and  the 

calculation  of  P^,  from  P^  i  =  2,  ...  p-1  requires  2^  ^2^2^*"^""^"  = 
computations.  Finally,  the  calculation  of  P  from  P  ^  requires  2^  2^  =  2^ 
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computations.  Thus,  the  total  number  of  computations  is  given  by 


^  22p+l-l  +  23p  =  1  ^(3^2) 

1= 1 


For  large  N,  this  figure  is  roughly  a  factor  of  ^  N  smaller  than  the  number 
of  computations  required  by  straightforward  evaluation  of  equation  (2). 

A  further  reduction  in  the  number  of  computations  can  be  effected  if  P 
need  not  be  evaluated  for  all  possible  values  of  its  arguments.  For  example, 
assume  that  P  is  to  be  evaluated  for  all  velocity  resolution  cells  but  only 
for  the  M  smallest  acceleration  cells  where  M  is  of  the  form  M  =  2^)+e>, 

0  <  g  <  p.  (The  reason  for  assuming  M  to  be  of  this  form  will  become  apparent 
in  a  moment. )  In  this  case,  the  binary  expansion  for  m  requires  only  p  +  g 
instead  of  2p  binary  digits;  i.e.  m  =  mQ  +  ...  +  m^^  ^  2^+{j~\  Examination 
of  equations  (4),  (5),  and  (6)  now  reveals  that  the  number  of  computations 
required  for  P^  is  equal  to  2‘^>+^  for  lei,  ...,  g  and  equal  to  2^>+G-^'  for 
l  =  g  +1,  ...  p.  It  follows  that  the  total  number  of  computations  C„  is 
given  by 


M 


22p+l-l  +  (p_g)  22p+g-l 


1=1 


=  N(M-N)  +  lo^  (  fp) 


(7) 
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It  is  interesting  to  compare  the  value  of  given  by  equation  (6)  with 
the  number  of  computations  required  by (two  other  methods)  for  evaluating  P 
for  N  velocity  resolutions  cells  and  M  acceleration  resolution  cells.  Straight¬ 
forward  evaluation  of  equation  (2)  requires  NM(N-l)  computations;  thus  the 
efficiency  of  the  above  proposed  method  can  be  assessed  by  evaluating  the 
ratio 


M 

NM(N-l) 


1  -  S 

M 

N-l 


1 

2 (N-l) 


10S2  (  !-) 


(8) 


As  a  numerical  example,  consider  the  numbers  N  =  M  =  32  for  which  equation  (8) 


M 


yields  =  0.08.  This  illustrates  the  considerable  computational 

advantage  the  proposed  method  has  over  straightforward  evaluation  of  equation 

(2). 


Another  way  of  calculating  P  for  N  velocity  resolution  cells  and  M 


acceleration  resolution  cells  is  to  combine  the  acceleration  factor 


with 


the  data  r,  in  equation  (2)  and  then  apply  the  Cooley- Tukey  method  for  a  pure 
velocity  periodogram  for  each  desired  value  of  m.  This  approach  results  in  a 
total  of  NM  log^N  computations  which  when  compared  with  yields 

^  ~  M  +  2  loS2  ^  M  ^ 


M 


NM 


loggN 


loggN 


(9) 


Substituting  N  =  M  =  32  in  equation  (8)  results  in  C^/NM  log^  N  =  l/2  which 
means  that,  in  this  case,  our  method  is  only  a  factor  of  two  more  efficient 
than  the  modified  Cooley-Tukey  method. 
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